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Generalization of the one-dimensional totally asymmetric exclusion process (TASEP) with open 
boundary conditions in which particles are allowed to jump I sites ahead with the probability 
pi ~ l/l a+1 is studied by Monte Carlo simulations and the domain- wall approach. For a > 1 
the standard TASEP phase diagram is recovered, but the density profiles near the transition lines 
display new features when 1 < a < 2. At the first-order transition line, the domain- wall is localized 
and phase separation is observed. In the maximum-current phase the profile has an algebraic decay 
with a cr-dependent exponent. Within the a < 1 regime, where the transitions are found to be absent, 
analytical results in the continuum mean-field approximation are derived in the limit a = — 1. 



I. INTRODUCTION 

Asymmetric exclusion process (ASEP) represents one 
of the basic models in studyingthe rich behavior of sys- 
tems held far from equilibrium 1] . It was first introduced 
m 0- El as a lattice model for protein synthesis and since 
then has been related to a wide range of other nonequilib- 
rium phenomena, like the surface growth 4, 5] belonging 
to the KPZ universality class [(| or traffic flow 0. In one 
dimension, subject to open boundary conditions, ASEP 
exhibits nonequilibrium phase transitions of the first and 
second order, which cannot occur in one-dimensional sys- 
tems that are in thermal equilibrium and include only 
short-range interactions. The finding of the exact ana- 
lytical solution to this model jl, 0, EJ contributed to a 
better understanding of general mechanisms leading to 
phase transitions in systems out of equilibrium pdj . A 
number of different extensions of this model have shown 
that the ASEP phase diagram is quite robust to various 
modifications |12| . 

The formulation of a large-deviation function [pH Ip4l | , 
as a nonequilibrium equivalent of the free energy, pointed 
out the importance of the effective long-range interac- 
tions appearing in these processes. In the present pa- 
per we include the long-range effect into this model in 
an explicit way, and examine the consequences for the 
phase diagram and critical behavior at the transition 
lines. We investigate the generalization of TASEP in 
which long-range interactions are introduced by letting 
particles overpass each other by jumping any distance 
1 < I < L with the probability pi ~ Z _CT_1 . Similar ex- 
tensions have been considered in the high-speed Nagel- 
Schreckenberg traffic model 

urn 

where cars are allowed 
to speed up if there is more than one empty site in front. 
However, these models do not allow cars to pass each 
other, which is justified for modeling single-lane traffic. 

An example of transport where particles are allowed 
to overpass each other may be found in the movement of 
the genome regulatory proteins along the DNA. Namely, 
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proteins bind to non-specific sites at DNA and search for 
a specific target site either by one-dimensional diffusion 
along the DNA (sliding) or by continuous dissociation 
from and re-association to DNA. Once dissociated, pro- 
teins undergo three-dimensional diffusion and eventually 
rcbind to another site. This process is either macroscopic 
(fully dissociated protein rebinds to a random site on 
DNA) or microscopic (dissociated protein remains within 
the range of the electrostatic potential of DNA and re- 
binds to one of the several nearby sites) |16| . In terms 
of TASEP, the first case corresponds to the creation and 
annihilation of particles in the bulk with constant rates 
(l7| . In the second case, dissociation and re-association 
events are correlated in space, which can be described 
in TASEP through the possibility for particles to travel 
larger distances in a single jump. 

The paper is organized as follows: in Section II we 
describe the generalization of ASEP to the long-range 
hopping, along with the definition of the open bound- 
ary conditions adjusted to include the hopping with ex- 
tended range. In Section III we present numerical re- 
sults obtained by Monte Carlo simulations and compare 
them to the expected phase diagram obtained from the 
current-density relation (fundamental diagram). In Sec- 
tion IV we proceed to investigate density profiles. Using 
the domain-wall approach, we determine the criterion for 
the phase separation in terms of the parameter a. In the 
maximum-current phase, using the appropriate scaling of 
density profiles, we find that the density profiles decay al- 
gebraically with the exponent that depends continuously 
on (j for 1 < a < 2, while the usual exponent ^ from 
ASEP is recovered for a > 2. Conclusions are given in 
Section V. 



II. LONG-RANGE HOPPING WITH 
EXCLUSION 

We consider the totally asymmetric exclusion process 
(TASEP) on a one-dimensional chain consisting of L sites 
where each site i is either occupied by a particle (ji = 1) 
or empty (r^ = 0). Instead of only the nearest-neighbor 
hopping, we allow each particle to jump from the site i 
to any empty site j > i on the right with the probability 
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pi which decays with distance I 



J\ as 



1 



1 



pi = 



Ct(<7+l)I" +1 



(1) 



where the normalization Ci(<? + 1) = Yli=i * _<T l is the 
partial sum of the Riemann zeta function. In the limit 
a — > oo, the jump probability Q becomes pi = Si t i 
and the model reduces to the standard TASEP with the 
nearest-neighbor hopping. 

In order to define open boundary conditions we first 
introduce the particle reservoirs on the left and right end 
of the chain, with constant densities a and 1 — f3, respec- 
tively, comprising L sites each. Since each particle can 
travel up to L sites in one jump, we allow any particle 
from the original chain to jump to the right reservoir. 
For example, the particle at the site L can jump to any 
site in the right reservoir, whereas the particle at the first 
site can jump only to the site L + l (first site of the right 
reservoir). Therefore, the total probability that the 
particle will leave the chain from the site i is the sum of 
probabilities of all possible jumps to the right reservoir 
from the site i, multiplied by the probability (3 that a site 
in the right reservoir is not occupied: 



state C — {t}, let C 4 '- 7 be the configuration obtained 
from C, with n and tj exchanging places, and C l the 
configuration obtained from C with n — > 1 — Tj. The 
system then evolves in time in two ways: either by ex- 
changing particles and holes with the reservoirs (C — > C l ) 
or by exchanging the particle at the site i with the hole 
at the site j > i (C — > C 1 ^). The impact of each of 
these two contributions depends on the interaction pa- 
rameter a. To each update mechanism we can assign a 
characteristic length defined as the first moment of the 
corresponding jump probability distribution. 

For the jumps within the chain, we get the average 
distance Xl as 



i=i 



a(f +i)' 



(•5) 



Similarly, the average distances 7_ and 7+ from the 
boundaries at which the particles are injected and re- 
moved from the chain, respectively, are given by 



7_ = anx(er), 
where 7l(<t) is defined as 



7+ = /?7l(ct), 



(6) 



(3 



L 

£ 

j=L-i+l 



er+1 1 



(2) 



Similarly, only the particles from the left reservoir that 
are within a distance L from the site i in the chain can 
jump to the site i. Since the probability to find the par- 
ticle at any site in the left reservoir is a, the total proba- 
bility ai that the particle from the left reservoir will enter 
the chain at the site i is given by 



Cl( 



L 1 
a \ 1 

' 3=i J 



(3) 



Due to the particle-hole symmetry that is present in our 
model, the probabilities on and Pi are not independent. 
Namely, particles jumping to the right can be seen as 
holes jumping to the left, which gives the relation ctj = 
a/?L-i+i//3. 

Within the continuous-time dynamics, the time evolu- 
tion of the system is described by the master equation: 
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1=1 



lai 



s5> 



Xl(c 



1=1 



1 



a(a-l) 



Cl(ct) 



By comparing these two lengths in the limit L — * 00, 
two different regions may be distinguished with respect 
to the parameter a: for a > 2, both A = liniL^oo Xl and 
7 = limi^oo 7l are finite, whereas for 1 < a < 2, A is 
finite, but 7 diverges. In the first case the particles, on 
average, enter and leave the chain within some finite dis- 
tance 7 from the boundaries. In the second case, where 
7 is infinite, particles are, on average, created and anni- 
hilated in the bulk and the model is in the regime similar 
to the TASEP with the bulk reservoir [T^]. The results 
presented below confirm this simple estimate. 

From the master equation J3J, the time change of the 
average particle density (r,-(f)) at the site i is given by 
the discrete lattice continuity equation, 



j:(n(t)) =Ji- J i+ i, 



(8) 



C)P(C, W{C -» C')P(C, t) 



c 



(4) 

where P(C, t) is the probability of the system being in the 
state C (given by the particular configuration of particles 
{n\i = 1, L}) at time t and W(C -> C')At {At -> 0) 
is the probability of the system going from the state C 
to the state C in the time interval [t, t + At}. Given any 



where Ji is the total current of particles, both those pass- 
ing over the site i and those leaving it 

L i L i 

Ji= Yl a k(i-(n))+J2 £ Pi-k{Tk(i-n))+520k{n) 

fe=i+l fe=l l=i+l fc=l 

(9) 

In the stationary state, the current is site- independent, 
and we can simply write the total current Ji = Jl(o!, /3, a) 
as a function of parameters a, f3, a, and L. 
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In order to make a comparison with TASEP, we first 
examine the dependence of the total current on param- 
eters a and (3, and possible discontinuities in its deriva- 
tives leading to noncquilibrium phase transitions. 



III. PHASE DIAGRAM 

From the way the boundary conditions were con- 
structed, we already know that the special choice of a 
and (3 such that a — 1 — (3 is equivalent to the case where 
the boundary conditions are periodic. Since the trans- 
lational invariance is then restored, all configurations C 
are equally probable 4] . The stationary density profile is 
flat, i.e. (ri) = a, which upon inserting into the expres- 
sion for the current © gives 

j L (a,f3,a) =X L (a)a(l-a), a + (3=1. (10) 

In the case of periodic boundary conditions (r; = Ti+i,), 
the number of particles TV is conserved, and the average 
density profile (r^) is constant and equal to N/L. This 
gives the current-density relation (Fig. ^| 

3 {p) = X LP {l- P ). (11) 
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FIG. 1: Fundamental diagram of TASEP with long-range 
hopping with a = 1.5 (solid line) and TASEP with the 
nearest-neighbor hopping (dashed line). 

It has been argued [l8L , that the knowledge of the 
current-density relation j(p) is sufficient to determine the 
phase diagram of a driven diffusive system coupled to the 
left and right reservoirs with constant densities p_ = a 
and p+ = 1 — /3, respectively. Explicitly, it was suggested 
|19|| that the stationary current of the system with open 



boundaries, in case of a positive drift (j > 0), is always 
selected according to the extremal principle 

( max{j(p),p + < p< p-} p-> p+ 

J = \ • (12) 

[ min{j(p),p- < P < P +} P - < p+. 

In the present model, the current-density relation 111|) 
has several important properties: (i) except for the fac- 
tor A l (c) , which accounts for the average length of the 
jumps, it is the same as in TASEP; (ii) in the infinite 
system the factor A, and therefore the current as well, 
are finite only for <j > 1; (iii) j(p) is symmetric to the 
particle-hole exchange. From the last property it fol- 
lows that the phase diagram is symmetric to the mutual 
exchange of the parameters a and (3. If the extremal 
principle applies to our model, the first property means 
that the structure of the phase diagram is the same as 
in TASEP. The second property then ensures that the 
current is still well defined for an infinite system. 

In the present model, the current has only one max- 
imum at p* = 1/2, like in TASEP. Consequently, the 
expected phase diagram has three phases: a maximum- 
current phase in the domain p_ > p* , p + < p* , with the 
bulk density p — p* , a low-density phase in the domain 
p_ < p* , with the bulk density p = and a high- 
density phase in the domain 1 — p + > p*, with the bulk 
density p = 1 — p + . 

The agreement of these predictions with the Monte 
Carlo simulations for different parameters a, (3, and a is 
presented in the following section. 



A. Numerical results 

Monte Carlo (MC) simulations were performed with 
the random sequential update. In other words, in each 
discrete time step (L such steps making one Monte Carlo 
Step per site (MCS/site)), the site i is randomly cho- 
sen from the chain of L sites. If the site i is empty, a 
particle jumps from the left reservoir to the site i with 
the probability on. If the site i is occupied, an integer 
m € {1, . . . , L} is drawn from the probability distribu- 
tion p m and two cases are possible: in the first case, if 
the site i+m < L is empty, the particle at the site i jumps 
to the site i + m, whereas in the second case, if i + m > L, 
the particle leaves the chain with the probability (3. In 
a typical run of the simulation we calculated the time 
averages of density profiles and the particle current j'l, 
omitting the results from the first to MCS/site, where to 
is the time required for the system to reach the (quasi) 
stationary state. In exploring the phase diagram, a mod- 
est size L = 200 was sufficient, and we used t ~ to ~ 10 6 
MCS/site. 

We found the same phase diagram as in TASEP 
only for a > 1, which is in agreement with the fact 
that A diverges for a < 1. Since these two regions are 
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essentially different, we treat them separately. 

The case a > 1. For a > 1 the results of the MC sim- 
ulations, depending on parameters a and 0, reveal three 
phases with the boundaries (Fig. |2J) and bulk densities 
of each phase the same as in TASEP 0, [^} . In order 
to check the order of the phase transitions across the 
boundaries, we consider the bulk density p = (N)/L and 
the current jl as functions of a for two characteristic 
constant values of (3 (Fig. |3 Fig. ^J. Since the phase 
diagram is symmetric to the exchange a <-> /3, this is suf- 
ficient to cover the phase diagram in the whole region of 
the parameters < a, (3 < 1. 
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FIG. 2: Phase diagram for a > 1 from the results of Monte 
Carlo simulations consisting of low-density (A) , high-density 
(B) and maximum- current (C) phases. 

In Fig. |31 we observe that for a > 1 the bulk density 
has a discontinuity at the coexistence line a = (3 
between the low-density and high-density phases, which 
is characteristic of a first-order transition. In Fig. 
0] the continuous transition to the maximum-current 
phase is presented. Although the finite-size effects in 
the maximum-current phase are visible, this diagram 
is in agreement with the extremal principle (|12|) for 
p* = 1/2 and symmetric to the mutual exchange of 
parameters a and (3, as concluded above from symmetry 
considerations. 

The case a < 1. In this domain, the current diverges in 
the limit L — > oo and it turns out that this completely de- 
stroys the phase diagram. The results from Monte Carlo 
simulations show the absence of the first- and second- 
order phase transitions, as shown in Fig. [3] (for varying 
a at constant [3 = 0.3) and Fig. 0] (for varying a at 
constant (3 = 0.7). 

In the limiting case a = — 1, a simple analytical so- 
lution of the stationary Eq. JSJ may be derived in the 
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FIG. 3: Dependence of the bulk density p on a for j3 — 0.3 
and various a (L = 200, t = 2 ■ 10 6 MCS/site) compared to 
the bulk density in TASEP (dashed lines). 
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FIG. 4: Dependence of current ji (a, (3, (t)/\l on a for f3 = 0.7 
and various a (L = 200, t = 2 ■ 10 6 MCS/site) compared to 
the current in TASEP (dashed line) . 



continuum mean-field approximation for a = (3 (see Ap- 
pendix A); i.e., by assuming (TjTj) = (Tj)(Tj) and setting 
x = i/L in the limit L — > oo. In particular, the current 
can be written in the closed form as 



j/L = 



a(2a + 3) 



1 - 



2a 



: Arsh> 



I 1 -2a 
2a(2a + 3) 



(13) 

It should be pointed out that all derivatives of (|13fl are 
smooth functions of < a = j3 < 1. 



IV. DENSITY PROFILES 

Monte Carlo simulations show that in the region a > 1 , 
different phases have the same bulk density as in stan- 
dard TASEP. On the other hand, at the phase boundaries 
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the spatial dependence of density profiles differs from the 
one in TASEP and in general depends on a. The exact 
solution of the present model is unknown and the mean- 
field approximation does not simplify it since it gives the 
complicated recursion equation for the average density 
(ri) that includes all the other densities (ry), j ^ i. The 
only case we managed to solve in closed form is the lim- 
iting case a = — 1 in the continuum limit (see Appendix 
A for a = /?), which is out of our main scope due to the 
absence of phase transitions in the region a < 1. 

In this section, we show how the domain-wall approach, 
previously used in TASEP [2l|, can be applied to the 
present model for the special choice of parameters a = 
f3 < 1/2 as long as a > 1 and treated analytically in order 
to compare it with the results from the simulations. 

In the maximum-current phase, to which this approach 
does not apply, the spatial dependence of the density 
profile is investigated by applying finite-size scaling to 
the results obtained from the MC simulations. 



A. Line a = f3 

The domain-wall approach 0| in TASEP is a way to 
describe the collective motion of particles with a step-like 
shock profile performing the homogeneous random walk 
on a one-dimensional lattice consisting of L + 1 sites. The 
approach is based on the fact that the stationary conti- 
nuity equation JSJ in the hydrodynamic limit, which is 
precisely the inviscid Burgers equation |20| , has a multi- 
valued solution at some point x s - This can be described 
by the formation of a discontinuity (shock) in a density 
profile moving with the velocity v s |22(, 



j[p( x f)} ~ j[p(x s )] 

p(xf) - p{Xs) 



(14) 



Depending on the sign of this velocity, the shock in 
TASEP is localized either at the left boundary (high- 
density phase) or at the right boundary (low-density 
phase). For a = f3 the velocity v s vanishes and since 
the shock has equal probability to be at any site, the 
resulting density profile is linear. 

In order to apply this approach to the present model, 
one has to assume a particular shock profile. In gen- 
eral, this can be done by investigating the continuum 
mean-field approximation of the continuity equation (jSJ) 
and constructing the shock profile from the two solu- 
tions, each matching one boundary condition. Instead, 
for a > 1 , we approximated the instantaneous shock pro- 
file p s (x) with a step function, like it has been done in 
TASEP: 



This is justified by the fact that the phase diagram, deter- 
mined by the bulk density in each phase, is the same as in 
TASEP. However, the feature in which the present case 
differs from TASEP is that the hopping rates of the ran- 
dom walker are site-dependent, because the excess cur- 
rent of all particles entering and leaving the low-density 
or the high-density domain depends on the position of the 
shock. This feature was already encountered in TASEP 
with the "Langmuir kinetics" [24J , in which particles are 
created and annihilated in the bulk with constant rates 
flc and Cljj, respectively. 

The hopping rate ri from the site i to the site i + 1 
is proportional to the total current of particles entering 
and leaving the high-density domain p + = 1 — (3: 



1 - a - (3 



£>(1 - /3) - 



i = 1, ...,L. 



(16) 



Similarly, the hopping rate U from the site i to the site i — 
1 is proportional to the total current of particles entering 
and leaving the low-density domain p_ = a: 



1 - a - 



i-1 

E 



(3 ) a 



i = 2,..., L+1. 

(17) 
h = and 



The shock is reflected at the boundaries, i.e. 
rL+i = 0. 

If we compare the bias Vi = r»— Z» of the random walker in 
the present model with the one from TASEP (vi = (3— a), 
we see that, as long as a 7^ /?, the qualitative picture from 
TASEP remains the same: for a < (3 we have Vi > 0, Vi 
and the random walker is localized at the right end (low- 
density phase); for a > (3 we have v.^ < 0, V« and the 
random walker is localized at the left end (high-density 
phase). However, on the coexistence line a = (3 the bias 
is not zero, except at the specific site % = L/2+1 {x « 1/2 
for L> 1), where both hopping rates and h are the 
same. 

The stationary probability Pi to find the random 
walker at the site i is given by the set of equations [23| 

n^Pi-i + h+iPi+i - (n + k)Pi =0, i = 2, L 

l 2 P 2 - nPi = 

r L P L + -Il+iPl+i = 0, 

(18) 

where the last two equations account for the reflecting 
boundary conditions. The solution to the above equa- 
tions l|18|) can be expressed in the closed form as 



a 



Ps(x) = 



1 — f3 x > x s 



(15) 



Pi 



Zl 



Z L 



(19) 
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FIG. 5: (Color online) Comparison of the density profiles (a = 
13 — 0.2, a — 1.5) for various chain sizes L = 10, 100, and 
1000 obtained from MC simulations (symbols) and from the 
domain wall approach (solid lines). 



FIG. 6: (Color online) Comparison of the density profiles (a = 
/3 = 0.2, a = 3.0) for various chain sizes L = 10, 100, and 1000 
obtained from simulations (symbols) and from the domain- 
wall approach (solid lines). 



where Zl is the normalization. In the last expression we 
have rewritten Pj in another form using the potential Vi , 



i-l 



L + l > i > 1. 



j'+i 



(20) 



Once the probability Pj is known, the density profile 
( T i) DW in the domain- wall approximation is simply given 
by M 



\DW 




(3). (21) 



The expression H21(l can be evaluated numerically and 
compared with the profiles obtained by MC simulations. 
In Fig. 0this was done for a = 1.5 and various system 
sizes (L — 10, 100, and 1000). The agreement between 
the MC data and Eq. <|21|l is excellent and justifies our 
initial assumption, even for very small system sizes (L 
10;. (A slight discrepancy between Eq. i|21[) and the 
actual profile obtained by MC simulations was observed 
only for a close to 1 , and it could be reduced by increasing 
the size of the system.) 

The nonlinear shape of the profile becomes more pro- 
nounced as L increases. For example, in Figure El the 
difference between the density profile for a = 1.5 and the 
linear profile is more evident for larger system sizes, sug- 
gesting that in the thermodynamic limit the shock tends 
to localize in the middle of the chain. 

As we increase the parameter a, the deviation from the 
linear profile is less visible. In the limit a — > oo, which 
is exactly TASEP, we should obtain the linear profile 
again. Already for a = 3.0 (Fig. HJJ) the density profile is 
almost linear for any system size, similarly as in the case 



of TASEP, and the aim is to find the threshold at which 
the standard TASEP reg ime sets on. 

As pointed out in |25| , the localization of the random 
walker is possible when the effective potential (|20f) is not 
homogeneous but has minima and maxima. If the poten- 
tial has only one minimum, as in the present model, the 
localization is possible in the sense that 



lim 

L — >oo 



[<z 2 ) s - (x) 



211/2 



= 0, 



(22) 



where x is the position of the random walker and (•••)« 
denotes the average taken over the stationary probability 
Pi . Using the above definition we plot the standard devi- 



ation Al = [(a 



J 1 / 2 as a function of the system 



size L in Fig. [7| We find in good agreement with the 
power law 



Li 1 < cr < 2 
L a > 2 



(23) 



Therefore, the domain wall is localized for 1 < a < 2 and 
delocalized for a > 2. This result can be obtained ex- 
plicitly by taking the continuum limit of the expressions 
(O, C3), CU), and (EOJ and checking the Taylor expan- 
sion of the potential 1)20(1 around the minimum located 
at the site i = L/2 + 1 (for even system size L). Since 
for a — [3 the symmetry relation lj = r^-i+2 holds, only 
the even powers of [x — 1/2) are nonvanishing: 



V{x) = V(l/2) + c 2 {x~l/2) 2 + c 4 (x-l/2) 4 + - 



(24) 



It can be shown (see Appendix B) that both C2 and C4 
are of the same order in L, 



c 2 ~ c 4 



-1-0 



(25) 
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FIG. 7: Standard deviation of the domain wall position ac- 
cording to the stationary probability Pi for a — (3 — 0.1 and 



various a. 



Therefore, if we keep only the quadratic term in the Tay- 
lor expansion, we obtain the Gaussian distribution and 
the variance proportional to 1/ y/c2 yields exactly the re- 
sult (|23)> . 

The phase separation has already been encountered 
in TASEP with "Langmuir kinetics" , first for particle 
creation and annihilation rates £Ia and f^u that depend 
on the system size as L _1 Ww a nd later for Qa,£Id ~ 
L~ a as long as 1 < a < 2 26] . In the first case the 
width of the shock was found to decay as L -1 / 2 , whereas 
in the second case the decay was still algebraic but with 
the exponent a/2. Since the similarity between these 
two models and the one we consider here lies in the fact 
that both models have the bulk reservoir, we checked the 
dependence of the rates oti and (3i on the system size L 
for er > 1 (see Appendix B). In the continuum limit and 
away from the boundaries, we found a similar power-law 
dependence on the system size, i.e. a(x),j3{x) ~ L~ a , 
so that the comparison can be made with TASEP with 
"Langmuir kinetics" and where a = a. In terms of the 
characteristic length 7 Q, it is clear that the similarity 
between the present model and the other two models with 
bulk reservoirs is established only when 7 diverges, since 
in this case the boundaries have a long-range influence 
on the bulk dynamics. 



B. Maximum-current phase 

As discussed in Sec. Ill (Cf. Fig. |2J, the maximum- 
current phase also appears in the presence of long-range 
hopping when a > 1 and sets in by a continuous phase 
transition for the same values of parameters (a, (3) as 
in the standard TASEP. However, one may expect that 
long-range hopping would alter the long-range correla- 
tions and therefore the behavior at the second-order 
phase transition to the maximum-current phase. 



Cl 
< 



0.2 
0.1 


-0.1 



-0.2, 




0.2 0.4 0.6 0.8 
k/L 




0.4 0.6 
k/L 



FIG. 8: (Color online) Deviation Ap(k/L) of the density pro- 
file from its bulk value p = 1/2 for a — 1.5 and various system 
sizes L (a — (3 — 0.8), before (a) and after (b) the scaling re- 
lation 12611 was applied. 



The exact solution in the case of TASEP with the 
nearest-neighbor hopping 0, 0| shows that the transi- 
tion to the maximum-current phase is characterized by an 
algebraic decay of the correlation function, while the cor- 
relation length remains infinite in the entire maximum- 
current phase. This is reflected in the algebraic decay of 
the density profile toward its bulk value p = 1/2 with the 
exponent 1/2. 

We studied this decay in the presence of long-range 
hopping by considering the deviation of a profile from its 
bulk value, Ap(fc) = |(rfc) — 1/2|, by using Monte Carlo 
simulations. We found that for given a, /3, and a and for 
different system sizes L\ and L 2 , Ap(fc) obeys the scaling 
relation 



Ap{kjL x ) = 



L2 
Li 



Ap(k/L 2 ), 



(26) 



where x depends continuously on a for 1 < a < 2 and 
is equal to the exponent 1/2 for a > 2. Fig. [S^ shows 
typical results of Monte Carlo simulations (a — 1.5, a = 
(3 = 0.8) for Ap(k/L) for various system sizes L. The 
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FIG. 9: (Color online) Deviation Ap(k/L) of the density 
profile from its bulk value p = 1/2 for various system sizes 
L = 800, 1500, and 5000 (a = @ = 0.8) for a = 1.2, 1.5, and 
1.8 (from top to bottom). The profiles Ap(k/L) for the same 
er are scaled to the profile Ap(k/Lo) (Lo = 5000) according 
to Eq. J2>J. 



but similarly to the model which includes Langmuir ki- 
netics m. 

Like in the standard TASEP, the second-order phase 
transition to the maximum-current phase has an infinite 
correlation length (as long as a > 1), which remains infi- 
nite in the whole maximum-current phase, and the profile 
displays an algebraic decay. However, the corresponding 
exponent is different from 1/2 in the region 1 < a < 2, 
where it depends on the range parameter a. 

In the special case a = — 1, which belongs to the regime 
without phase transitions, a < 1, and corresponds to uni- 
form hopping probability, we derived the exact analytical 
results both for the bulk properties and for the density 
profile along the line a — (3. 

The dynamical scaling, now under consideration, 
should give a more complete description of the consid- 
ered second-order phase transition. It would also be of 
interest to examine the relation to some other models for 
nonequilibrium phase transitions with long-range inter- 
actions, in particular the models including temperature 
such as the two-temperature kinetic Ising model. 



best numerical fit to the scaling Eq. i|26|) . for which all 
data collapse on a single curve (see Fig. HJd), is achieved 
for x = 1/4. In Fig. [5] we present the collapsing fits 
for several values of a. The best numerical results for 
collapsing fits are achieved for x = (er— 1)/2 which, along 
with x = 1/2 for a > 2, leads to the conjecture 



x = min{— — ,-), <j>\. (27) 



V. CONCLUSION 



APPENDIX A: CASE a = -1 IN THE 
CONTINUUM MEAN-FIELD APPROXIMATION 

In the continuum limit we set i/L — > x and substitute 
all sums in Eq. © with integrals. When er = — 1, we get 
simple expressions for the probabilities a(x) and (3(x): 

a(x) = a(l-x), j3(x) = fix. (Al) 

Inserting (|A1|) in (JSJ) gives the first-order nonlinear dif- 
ferential equation 



We considered the totally asymmetric exclusion pro- 
cess (TASEP) with open boundary conditions, general- 
ized in such a way that the long-range effect was intro- 
duced directly, through long-range hopping, where the 
probability of distant hoppings decays as a power law 
Boundary conditions were adjusted in such a 
way that the periodic boundary conditions are attained 
in the special case a = 1 — /3. 

We found that for a > 1 the phase diagram of this 
generalized model remains the same as the one in TASEP 
with the nearest-neighbor hopping, containing the same 
first- and second-order transition lines. The bulk density 
equals the one in the standard TASEP, while the current 
is equal to the one in TASEP normalized by a factor- 
depending on er. 

As expected, new effects were found at the transition 
lines, where long range of hopping has an impact on long- 
range correlations which determine the character of the 
transitions. They were observed by studying the density 
profile. 

At the first-order transition line a = (3 the profile be- 
comes localized, in contrast to the standard TASEP case, 



dn 
dx 



[2n+((3-a-l)x+(l-n+a)]-n-a(l-x) = (A2) 



with the boundary conditions 



n(0) = and n(l) = /i. 



(A3) 



In the above expressions, 4| = p(x), and pi(a,(3) = 

Jo p(x)dx is the bulk density. For a = (3 the first-order 
differential equation (|A2|) is exact, leading to 



— |[n(a;)] 2 + n(x) (— + a — x^j + ^—ax 2 — ax^j X = 0, 

(A4) 

where we used the fact that fi(a, a) = 1/2. From (|A4() 
the density profile is easily obtained as 



p(x) 



1 1 



(1 - 2a) (as - 1/2) 



2 2 ^2(1 _ 2a ) - x (l - 2a) + (a + 1/2) 2 ' 

(A5) 
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Finally, the total current j/L of particles entering or leav- 
ing the chain can be obtained by inserting l|A5[) into the 
expression for current, so that 

j/L = a [ (1 - x)[l - p{x)]dx = j3 [ xp(x), (A6) 
Jo Jo 

which gives Eq. 



APPENDIX B: ASYMPTOTIC EXPANSIONS OF 
an, Pi, n AND k IN THE LIMIT OF LARGE L 



1 — a — p 



f(l-x + 2e,a,(3, a) 



-,a-X 



(B6) 



where e = 1/L, and the function f(x, a, a, (3) is finite for 

Q < x < 1. We define A-™^ as the discrete n-th order 
"derivation" of Vi 



We start by rewriting © and (01 in a more convenient 
form: 



A (n) = A ( ?7 1} - A ( "~ 1} 



(B7) 



a; = a I 1 



Ci-i(g + i) 

a(^+i) 



A = j9 1 



a(a + i) 

(Bl) 

Assuming the convergence of the Zeta series, we estimate 
the partial sum in the limit of large L, 



a(<r+i) = c(a + i) 



oL c 



+ o(L-°- 1 ). 



(B2) 



After inserting i|B2|) in the expressions J3J and J2J and 
taking the continuum limit with x = i/L, we get for 
< a; < 1 



with Aj = K — Vi—\. In the continuum limit with 
x = i/L we get 



A 



(n) 



(B8) 



Finally, inserting the expressions l|16|) and (|17|l in the 
potential ifflty and expanding the function f(x,a,a, (3) 
around the minimum of the potential, x = 1/2 + e/2, we 
obtain the result for the first two expansion coefficients: 



^^^[(*-i/z)--i] 



C(^ [(1 -^- l] 



(B3) 



(B4) 



Similarly, the expressions for the hopping rates r- and U 
in the continuum limit for 0<i<1 are 



r(x) 



~ 0) 



f(x,a,a,/3) 



(B5) 



=1/2 



2f(l/2,q, q,q) £2 _^ 
Aa(l — a) 



(B9) 



C4 



2/'"(l/2,a,Q,o)_ i2 _„ 



X=l/2 



Aa(l — a) 



(BIO) 
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